Kinematic Segregation of Flowing Grains in Sandpiles 
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We study the segregation of granular mixtures in two-dimensional silos using a set of coupled 
equations for surface flows of grains. We study the thick flow regime, where the grains are segregated 
in the rolling phase. We incorporate this dynamical segregation process, called kinematic sieving, 
free-surface segregation or percolation, into the theoretical formalism and calculate the profiles of 
the rolling species and the concentration of grains in the bulk in the steady state. Our solution 
shows the segregation of the mixture with the large grains being found at the bottom of the pile in 
qualitative agreement with experiments. 



I. INTRODUCTION 

When a mixture of grains | differing in size and 
shape is poured in a two-dimensional cell consisting of 
two parallel vertical slabs separated by a narrow gap (i.e., 
a granular Hclc-Shaw cell) the grains segregate accord- 
ing to their size and shape. The typical experimental set 
consists of a vertical "quasi-two-dimensional" Hele-Shaw 
cell with a narrow gap of 5 mm separating two transpar- 
ent plates of 300 mm by 200 mm (see Fig. Q). We close 
the left edge of the cell leaving the right edge free, and 
we pour continuously, near the left edge, an equal- volume 
mixture of grains differing in size and shape. When the 
larger grains are more rounded than the small grains, 
a segregation of the mixture is observed with the large 
grains being found at the bottom of the cell and the 
small at the top. When the large grains are rougher than 
the small grains a periodic pattern consisting of layers of 
large and small grains parallel to the pile surface occurs 
|7H13| . This periodic pattern is called spontaneous strat- 
ification and it was found to occur for a wide range of 
size ratios between the grains: &ijd\ > 1.4, where d\ is 
the typical size of the small grains and di is the typical 
size of the large grains. When the size ratio between the 
grains is close to one, only segregation is found no matter 
the shape of the grains, i.e. no matter the angle of repose 
of the species. 

An important segregation mechanism acting when 
c?2 / d\ is not close to one, and in the the case of thick flows 
is the so called kinematic sieving! l3, free surface segre- 
gation [|l5| or percolation effect Due to this phe- 
nomenon the large grains in the rolling phase are found to 
rise to the top of the rolling phase while the small grains 
sink downward through the gaps left by the motion of 
larger grains in the rolling phase (Fig. |l|) . Since only the 
small grains interact with the surface they are captured 
at the surface first causing the larger grains to be con- 
verted further to the bottom of the pile. Thus, a strong 



segregation effect results in a pile with the small grains 
being found at the top of the pile and the large grains be- 
ing found at the bottom of the pile. This effect has been 
observed in high-speed photography experiments to be 
important for segregation and stratification of granular 
mixtures |s| in well developed flows of granular mixtures. 

In this paper, we study analytically the segregation of 
granular mixtures poured in two-dimensional sandpiles 
due percolation in the rolling phase by using a suitable 
modification to the equations of motion for the rolling 
species interacting with the static grains of the pile. We 
calculate the corrections to the steady-state profiles due 
to the percolation effect and compare with theoretical 
predictions without percolation effects and with experi- 
mental observations. 

Steady state solutions for the granular flows of grains 
in two-dimensional geometries have been calculated in 
]l(| for the case of a single species flows in thin rotat- 
ing drums, in |l7| for the case of two-species differing in 
angle of repose and flowing in two-dimensional silos, and 
in [p|[ll|l for thin flows of grains in two-dimensional si- 
los with two species of grains differing in size and shape. 
When the grains differ only in angle of repose, the solu- 
tion shows the segregation of the mixture and the con- 
centrations decay slowly as a power-law of the position in 
the cell 17 . When there is a large difference in size and 



shape the solution shows a strong segregation of the mix- 
ture and the profiles decay exponentially fast |] jll[] . The 
segregation solution is only valid for mixtures of large- 
rounded grains and small rough grains. Otherwise, the 
steady state solution is unstable, leading to the stratifi- 
cation of the mixture. 

Here, we treat the case when the difference in size of 
the species is not too small (i.e. d 2 /di > 1.4) but, as op- 
posed to we focus on the thick flow regime. In this 
regime the thickness of the rolling phase is expected to be 
large, and it is possible to observe the percolation effect 
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(the segregation in the rolling phase) ||. The steady- 
state solution we find shows the complete segregation of 
the mixture with the large grains being found at the bot- 
tom of the pile, and we compare this solution with the 
solution corresponding to the case of thin flows when the 
percolation effect is absent. 

The paper is organized as follows: In Section || we 
present the theoretical form alis m for surface flows of 
granular mixtures. In Section III we calculate the steady- 
state solution of the problem, discuss our results, and 
compare to the theoretical predictions obtained without 
the percolation effect. 



incoming flux 
y of grains 



percolation effect 
^of$ m tne rolling phase 




FIG. 1. Schematic representation of the percolation effect 
in the rolling phase. We pour continuously a mixture of large 
and small grain in a two-dimensional cell and a steady thick 
flow of grains is observed. The percolation effect consists on 
the size segregation of grains in the rolling phase: large grains 
are observed to rise to the top of the rolling phase, and small 
grains drifts to the bottom of the rolling phase. Due to the 
percolation effect, the small grains are the first to be captured 
at the pile surface resulting in the segregation of the mixture 
in the bulk. 



II. SURFACE FLOW OF GRANULAR MIXTURES 

The theoretical study of surface flows of granular ma- 
terials was triggered by the works of Bouchaud, Cates, 
R. Prakash, and Edwards (BCRE) (Ts| and Mehta and 
collaborators In a recent theoretical study for the 

case of a single-species sandpile BCRE Q proposed 
two coupled variables to describe the dynamics of two- 
dimensional sandpile surfaces: the local angle of the 
sandpile 9(x,t) (or alternatively the height of the sand- 



pile h(x,t)) which describes the static phase (i.e., the 
grains which belong to the pile), and the local thickness 
of the layer of rolling grains R(x, t) to describe the rolling 
phase (i.e., the grains that are not part of the pile but 
roll downwards on top of the static phase). BCRE also 
proposed a set of convective-diffusion equation for the 
rolling grains, which was later simplified by de Genncs 
II- 

Recently, Boutreux and de Gennes (BdG) |l7| have ex- 
tended the BCRE formalism to the case of two species. 
This formalism considers the two local "equivalent thick- 
nesses" of the species in the rolling phase R a (x,t) (i.e. 
the total thickness of the rolling phase multiplied by the 
local volume fraction of the a grains in the rolling phase 
at position x), with a = 1,2 respectively for small and 
large grains. The total thickness of the rolling phase is 
defined as 



R(x,t) = Ri(x,t) + R 2 (x,t). 



(1) 



The static phase is described by the height of the sandpile 
h(x,t), and the volume fraction of static grains <fi a (x,t) 
of type a at the surface of the pile. We consider a silo or 
cell of lateral size L and the pouring point is assumed to 
be at x = 0. For notational convenience we do not con- 
sider the difference between the angle and the tangent of 
the angle, i.e. 0(x,t) = — Jj|. The equations of motion 
for the rolling species are |l7f 



dR a (x, t) 



dR a , r 



dt a dx 

and the equation for h(x, t) follows by conservation 



(2a) 



dh(x,t) 
dt 



= -r x -r 2 . 



(2b) 



Here v a is the downhill convection velocity of species 
a along x, assumed to be constant in space and in time. 
The interaction term T a takes into account the conver- 
sion of static grains into rolling grains and vice versa, and 
it is defined as 



Ti^aiiO^xRi-bx^Ri 



T 2 = a 2 {9)4> 2 R2e- XRllR - b 2 {9)R 2 e 



-Afli/ii 



(3a) 



(3b) 



This definition involves a set of a priori unknown col- 
lision functions contributing to the rate processes: a a (9) 
is the contribution due to an amplification process, (i.e., 
when a static grain of type a is converted into a rolling 
grain due to a collision by a rolling grain of type a), 
b a (9) is the contribution due to capture of a rolling grain 
of type a, (i.e., when a rolling grain of type a is converted 
into a static grain). Cross-amplification processes, (i.e., 
the amplification of a static grain of type /3 due to a col- 
lision by a rolling grain of type a) do not contribute to 
the exchange of grains since these kind of processes are 
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inhibited by the percolation effect. Therefore, we do not 
consider them in the definition of the interaction term. 

The interaction between rolling and static grains is as- 
sumed to be proportional to the number of interacting 
rolling grain, then it is proportional to R a (x,t). This 
is clearly valid for thin flows. However, in the case of 
thick flows treated here, this approximation might still 
be valid [ [To[ since the interaction might be proportional 
to the pressure exerted by the fluid phase which in turn 
is proportional to thickness of the rolling phase |2(J • 

Due to percolation effect the interaction with the bulk 
of the large rolling grains R%(x,t) is greatly reduced as 
long as there are small rolling grain R\{x,t) interacting 
with the surface. Therefore, to take into account the 
percolation effect, we replace R2(x,t) in the definition 
of the interaction term by i?2(x, f) exp[Ai?i (a;, t)/R(x, t)]. 
The exponential factor multiplying R2(x,t) mimics the 
fact that the interaction of large grains is screened by 
the presence of small grains so that R2 interact with 
the grain at the surface of the static phase only when 
Ri(x,t) <C R(x,t)/X. The dimensionless parameter 
A > measures the degree of percolation, with A = 
corresponding to the no percolation case treated in Q . 

The concentrations of static grains at the surface of 
the pile 4> a (x,t) are given by 



4>a(x,t)^ = -F a , 



and 



>l -t- <P2 



= 1. 



(4a) 



(4b) 



We use the following definitions of the collision func- 
tions 101: 



a a (0) = l aa Il[6(x,t) - 6 a (<j>p)] 
b a {9) = laa n[6 a (<f> )-e{x,t)} 



where 



U\x] 



if x < 
x if x > 



(5a) 



(5b) 



Here, the rates ~f aa > has dimension of inverse time, 
and v a /j aa ~ d a . The generalized angle of repose 9 a ((f>p) 
of a a type of rolling grain is a continuous function of the 
composition of the surface 4>p |0J, and we also define 9 a p 
as 9 a ((f>p) for 0/3 = 1 



92) 



m4>2 
m4> 2 



in 

hi = -m(j)i 



'22, 



(6) 



where m = 9\2 — On — ^22 — #21- We assume the dif- 
ference tp = 9\{cj)2) — 6*2(02) to be independent of the 
concentration (f>2- We notice that for mixtures of grains 
with different shapes we have #11 7^ 6*22 (here 6>n is the 
angle of the pure small species and #22 denote the angle 
of repose of the pure large species), and different sizes 
6*12 7^ 6*21 ■ hi the following we treat only the case when 



the type 1 small grains are the roughest (i.e. 6*22 < #11, 
since the angle of repose depends on the surface proper- 
ties of the grains and it is larger for the roughest grains) 
or when the grains differ only in size (i.e., 6*22 = On). The 
case when the large grains are the roughest (6*11 < #22) 
leads to stratification of the mixture and it is treated in 

H3- 




e 2 (4> 2 ) 



FIG. 2. Plot of capture b a (9) and amplification a a (9) func- 
tions. These functions are expected to be continuous in a re- 
gion near the angle of repose (as shown by the thin curves). 
However, when the the ratio between the size of the grains is 
dild\ > 1.4 then tp = 61(412) — #2(^2) is large enough so we 
can approximate these functions by the forms shown in this 
figure. 

The collision functions a a and b a are expected to be 
continuous in the region of interest near the angle of re- 
pose (as shown, for instance, by the thin curves in Fig. |^). 
However, when the ratio between the size of the grains is 
large enough, the angular difference ip = 6*i(0 2 ) — 9 2 {(j)2) 
is expected to be large, and in this case, we can approx- 
imate the functions by the forms defined in (j^). When 
d 2 /di is not close to one (i.e. c^/^i > 1-4 according 
to experiments) strong segregation effects act in the sys- 
tem. In the case of thin flows treated in |)| (i.e. for 

R(x,t) ~ 3 * g?2) there is no percolation in the rolling 
phase and all grains (large and small) interact with the 
surface, but the large grains are not being captured be- 
cause of the large difference in size. Thus the small grains 
are the only ones to be effectively interacting with the 
sandpile surface resulting in a strong segregation effect. If 

the flux of grains is large enough (i.e., for R(x, t) ~ 3*c?2, 
the case treated here), segregation occurs in the rolling 
phase due to the percolation effect so that the large 
grains migrate to the top of the rolling phase and the 
small grains migrate to the bottom of the rolling phase 
|]i^ , ^5| ]8|JlCt| . Then small grains are the only ones inter- 
acting with the sandpile surface and again strong seg- 
regation is expected. This is confirmed by experiments 
done for ^2/^1 > 1-4 showing sharp segregation patterns 
when the large grains are smoother than the small grains 
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S = 822 — On < 0, and stratification when 5 > 0. 

When the size ratio is close to one (cfe/^i < 1<4), the 
angle ip is expected to be small and the collision func- 
tions a a and b a can be linearized around the respective 
angles of repose 8 a (4>2), since the region of interest is lo- 
calized to a small region around the angles of repose of 
the species. In this case segregation effects are expected 
to be less important in comparison with the above case. 
Indeed, according to the experiments (l^] mixtures of 
grains differing in size and shape with size ratio close 
to one {da/d-y < 1.4) give rise only to segregation and 
the concentration show slow decay, i.e. the segregation 
pattern is not sharp as in the above case. The case of 
weak segregation is treated theoretically in , where it 
is shown that the concentration profiles show slow (al- 
gebraic or power-law type) decay as a function of the 
position in the cell. 



III. PROFILES IN THE STEADY-STATE 

We now calculate the steady-state solution of the equa- 
tions of motion for the two-species sandpile including the 
effect of percolation. 

We consider the geometry of a silo of lateral size L. 
We assume that the difference if) = 81(412) — 82(42) is in- 
dependent of the concentration 42, then ip = 6\\ — 821 = 
812 — 822- We set v\ — V2 = v, and 711 = 722 = 7- We 
seek a solution where the profiles of the sandpile and of 
the rolling grains do not change in time. Since stratifi- 
cation is an oscillatory solution, stratification cannot be 
observed for the steady-state solution. We set 



and 



dh 

at 



8R a (x) 
dt 



R° 



= 0, 



with the following boundary conditions: 

R a (x = 0) = Rl, R a (x = L) = 0, 



(7a) 



(7b) 



(8) 



with R° = r% + r%. 

The profile of the total rolling species R(x) 
R2(x) is 

R° 

R(x) = ^-(L-x). 



Ri(x) 



(9) 



which is obtained from J2a|), (2b), and ((?]). The linear 
profile for the total number of rolling species was also 
found in the case of one type of species in Jig] , In fact 
this result is independent of the type of interaction term 
used, and it is a consequence of the conservation of grains. 

Th e eq uations for the rolling species are obtained from 
( pa| ) , d4a| ) , and (0) as a function of the concentrations 4a 



dR a (x) 
dx 



= -(vR°/L) 4> a (x), 



(10) 



Since the collisions functions are defined according to 
the value of the angle in comparison to the generalized 
angles of repose, we divide the calculations in two re- 
gions: Region A, where 82(42) < 8 < 81(42), and Region 
B, where 8 < 82(42) < 0i(<fe). 

Region A. If 8 2 < 8 < 6\, we obtain 



4i(x) 



1, 



42(x) =0. 



(11) 



then we obtain the profile of the rolling species using ( fTo| ) 
and (H) 

R x (x)=R\-? I -x, R 2 (x) = R%. (12) 
The profile of the sandpile is obtained from (|t]) , (||) , and 



8(x) 



-v/7 



LR Q 1 /R°-x' 



(13) 



This solution is valid when 8(x) > 81(42 = 1) = #21, or 
for x < x mi where x m is 



R° 



L 



lip' 



(14) 



We notice that in this region the solution is the same 
as the one found for the case of thin flows without per- 
colation (A = 0) in |f[ll| . This is because, even though 
the large grains may interact with the pile surface when 
percolation effects are absent, the small grains are cap- 
tured before the large grains since the surface at the top 
of the pile is made of small grains, so that the large grains 
roll down very easily and they are not trapped (in other 
words, the angle of repose of the pure small grains 6\i is 
larger than the angle of repose of the large grains rolling 
on a surface of small grains #21 )■ The profile of the small 
rolling grains decays linearly with the distance as in the 
case of the total rolling species Eq. (^), since essentially 
the small grains are the only ones effectively interacting 
with the surface, so that it is a single species problem. 
Next, we show that the percolation effect shows up in 
the region at the bottom of the pile causing corrections 
to the exponential profiles of the concentrations found in 

Si- 
Region B. If 9 < 8 2 < 9t, from ©, (§), and @ we 
obtain 



vR° 



= [8(x) - ex(42)]Rx[x) + [8(x) - 82(42)]R2(x)e- XR '/ R , 



and therefore 
8(x) - 1 



2[<P2) 



-vR°/(jL) +-4R ± (x) 
R*(x) : 



(15) 



(16) 
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where R*(x) = R 1 (x) + R 2 {x) e - XR ^ R . Using @ we ob- 
tain the concentration as a function of the rolling species 



(17) 



We obtain the equations for the rolling species using Eqs. 
© and (0), 



dR^x) _ gi(g) /g° flaCarje-^ 1 /-" 
9a; ~~ ~i?*(x) + r 



(18) 



where r = v/(^ip). We set 



m(x) 



i?i(x) 



i?*(x)' 
and we obtain from ( |lS| ) 

1 (1 - u)e" Al 



(19) 



u'(x) = —u I — 



dx R*(x) 



dR*(x) 1 
<9x i?*(x) 

(20) 



This equation cannot be solved in closed form, but to 
a good approximation we assume that dR*(x)/dx ~ 
dR(x)/dx (valid for small A), and the equation is 



u{x) w -u(l - u)e _Au /r, 



which has the solution 



Ri(x) 



R*(x) 1 + Ce-( x - x ^/ r 



(21) 



(22) 



where C is an integration constant obtained by consider- 
ing the continuity at X = x m . When there is no percola- 
tion effect (A = 0) we recuperate the result found in [ p"l| , 
i.e. Rjjx)/ R(x) ~ exp[— (x — x m )/r] when r C I. Using 
Eq. (XljY) the concentration profiles can be obtained as a 
function of position. 

To compare this solution with the case of no percola- 
tion effect, we obtain the exact numerical solution of Eq. 
( |l8| ) . Figure || shows the solution we find for the rolling 
species with percolation in the rolling phase in compar- 
ison with the solution corresponding to no percolation 
effect A = 0. We see that the percolation effect tends to 
increase the degree of segregation of the mixture in com- 
parison with the results found in |j) without percolation 
as expected; the decay of the small rolling species when 
there is segregation in the rolling phase is faster than the 
exponential decay found when there is no segregation in 
the rolling phase. 




percolation effect 
no percolation 



A(x) 



10 20 30 40 
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FIG. 3. Results of the profile of rolling species when the 
percolation effect is present. Notice that small species (type 
1) are the only one interacting with the pile up to the center 
of the pile. Then the large grains start to be captured to, 
and the number of large rolling grains R,2(x,t) decreases as 
well. For comparison we also plot in dashed line the solution 
corresponding to the case of thin flow where the percolation 
effect is not acting. We see that the segregation in the case 
of percolation is stronger than the segregation in the case of 
no percolation effect. 

The concentration of small grains decays very fast 
near the center of the pile as observed in experiments 
f|p"2"[ . with a region of mixing with a characteristic size 
of r = v/jip. For typical experimental values correspond- 
ing to a system of L = 30 cm, di = 0.27 mm, e?2 = 0.8 
mm we have B v ~ 10 cm/sec, 7 ~ 20/sec, R° ~ 0.5 
cm, tan?/; = tan#n — tan 6*21 — 0.2 [^2|. Therefore, the 
region of mixing in the center of the pile is of the order 
of u/(7tan-0) ~ 2.5 cm. 

So far we have calculated the steady state solution of 
the problem. Numerical simulations p- |TT| , ^l| shows that 
this type of solution is only stable when the small grains 
are rougher than the large grains. When the large grains 
are the rougher, the solution turns unstable and the sys- 
tem displays the stratification pattern found in . 



IV. DISCUSSION 

In summary, we study analytically segregation in gran- 
ular mixtures including the effect of segregation in the 
rolling phase, by introducing a suitable modification of 
the interaction term in the BdG equations. We find the 
corrections to the profiles of the steady-state of the filling 
of a silo. The profiles show strong segregation effects in 
agreement with experimental findings. The small correc- 
tions found here for the case of thick flows of grains with 
percolation in the rolling phase compared to the case of 
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thin flows but with a large difference in size ratio treated 
i n fll| suggests the applicability of the results found in 
pd| also to the case of thick flows. 
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